Quantum effects in a half-polarized pyrochlore antiferromagnet 
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We study quantum effects in a spin-3/2 antiferromagnet on the pyrochlore lattice in an external 
magnetic field, focusing on the vicinity of a plateau in the magnetization at half the saturation 
value, observed in CdC^CU and HgC^CU. Our theory, based on quantum fluctuations, predicts the 
existence of a symmetry-broken state on the plateau, even with only nearest-neighbor microscopic 
exchange. This symmetry broken state consists of a particular arrangement of spins polarized 
parallel and antiparallel to the field in a 3:1 ratio on each tetrahedron. It quadruples the lattice unit 
cell, and reduces the space group from Fd3m to P4332. We also predict that for fields just above 
the plateau, the low temperature phase has transverse spin order, describable as a Bose-Einstein 
condensate of magnons. Other comparisons to and suggestions for experiments are discussed. 

PACS numbers: 75.10.-b,75.10.Jm,75.25.+z 
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Frustrated quantum magnets, in which many 
symmetry-unrelated states are classically degenerate, 
are a fascinating venue in which to observe emergent phe- 
nomena. In most frustrated materials, the degeneracy is 
lifted classically, by lattice distortions (Jahn- Teller, spin- 
PeierlsQ), or longer-range (e.g. dipolar) interactions 0. 
More intriguing is the possibility that the degeneracy 
can be removed by quantum fluctuations. In such a 
"quantum order by disorder" scenario, the system would 
pick a ground state by maximally delocalizing among 
many degenerate classical states, thereby minimizing 
its quantum zero point energy. In toy models, this 
can lead to exotic ordered states (such as valence bond 
solids), or even "quantum spin liquid" states, where 
the classical degeneracy is split without any symmetry 
breaking. Experimentally, clean signatures of the lifting 
of degeneracy by quantum fluctuations have, however, 
remained elusive, presumably because of the dominance 
of the classical mechanisms discussed above. 

In this paper we describe a theoretical study of quan- 
tum effects in the insulating chromium oxides, CdCr2 04 
and HgCr204, in which magnetic spin-3/2 Cr 3+ ions form 
a pyrochlore lattice (a network of corner-sharing tetra- 
hedra). Because of the half- filled ti 9 magnetic levels, 
these materials lack orbital degeneracy, and are also very 
isotropic magnetically. An appropriate minimal theoret- 
ical model is thus the nearest-neighbor pyrochlore lattice 
Heisenberg antiferromagnet, with Hamiltonian 
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where Si is a spin-3/2 operator on site i of the pyrochlore 
lattice, {ij) indicates the sum is over nearest-neighbor 
bonds, and we have taken the field H to point along 
the z-axis. Most intriguingly, these materials display a 
robust low-temperature magnetization plateau in an ap- 
plied magnetic field, with magnetization well-quantized 
at half its saturation valueQ. While there has been a 




FIG. 1: (Color online) Planar projection of the "R" state, 
stabilized by quantum fluctuations. Dots indicate minority 
spins. The dashed lines are parts of the next layer of tetrahe- 
dra 



great deal of prior theoretical work on quantum effects 
m zero mag netic field @,@|, the connection to experiment 
seems unclear. Here we will instead focus on the physics 
on and near the plateau. First, we predict a magnetic 
structure on the plateau shown in Fig. with space 
group P4 3 32, which we call the "R" state. Second, we 
predict the existence of XY order (i.e. transverse to the 
field) at low temperature for fields just above the plateau. 
If quantum fluctuations dominate over second neighbor 
exchange, the transverse order is ferromagnetic. 

These conclusions obtain both for Eq. Q and a mod- 
ified model including spin-lattice coupling to local Ein- 
stein bond phonons, which Penc et ai.[jj have suggested 
play an important role. While we agree such coupling 
likely supports the large width of the plateau, our con- 
clusions with respect to the ordering symmetry are essen- 
tially independent of it. For clarity of presentation, we 
therefore suppress spin-lattice coupling in the following. 
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Our approach relies on the substantial polarization in 
the vicinity of the plateau, which indicates that, on aver- 
age, the transverse components of the spins are reduced. 
Moreover, the existence of a plateau with zero differen- 
tial susceptibility \ = dm z /dh = at T = implies the 
absence of transverse ordering, (Sf) — 0. We therefore 
split the Hamiltonian into terms involving longitudinal 
and transverse spin components, and treat the latter per- 
turbatively (though we will work to rather high orders). 
Specifically, H' = Ho +Hi, with 

n K s ? h f h2 ] ~ J * E (^) 2 . ( 2 ) 
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where Sf = Sf ± iSf are the usual ladder operators, 
Sf = Y)jcf Sj is the sum of longitudinal spin components 
on a tetrahedron labeled by t, and we have introduced the 
dimensionless magnetic field h = H/2J (we take h > 
without loss of generality). For bookkeeping purposes, we 
have written distinct longitudinal and tranverse exchange 
couplings J z , J± , though we are ultimately interested in 
J z = J± = J. We note that "easy axis" perturbation 
theory in a = J±/J z becomes increasingly accurate at 
larger h, and that the symmetry of H! is preserved for any 
a > 0. This gives us further confidence in the qualitative 
predictions of easy-axis perturbation theory. 

We now sketch our derivation of an effective model de- 
scribing the quantum mechanical lifting of the classical 
degeneracy. Consider first the th order Hamiltonian, 
Eq. (J2J. It is straightforward to see that Ho has three 
types of ground states, depending upon h: for h < 1.5, 
each tetrahedron has two spins with Sf — +| and two 
spins with Sf — — §; for 1.5 < h < 4.5, the half- 
magnetized region, each tetrahedron has 3 "majority" 
spins with Sf = +| and one "minority" spin Sf = — |; 
and all spins are fully polarized for h > 4.5. Therefore, 
this simplistic view predicts 3 magnetization plateaus ex- 
ist - zero, half and full magnetization. 

Focusing first on the plateau of interest, the th order 
"3:1" ground states are extensively degenerate, due to the 
freedom to locate the minority spin of each tetrahedron. 
We may view the transverse terms in Eq. as inducing 
"quantum fluctuations" into these states; in particular, 
the action of Hi transfers a total spin AS Z — ±1 between 
nearest neighbor sites. Technically, this must be analyzed 
by Degenerate Perturbation Theory (DPT). 

We define the projection operator, V, onto the ground 
state plateau subspace. Consider any exact eigenstate 
I*). Its projection |\J/ ) = V\^l) satisfies the "effective 
Schrodinger equation" 



VHx 



where the resolvent operator Q = e } Hq (1 — V) Hi is lin- 
ear in J±. Because the resolvent contains the exact en- 
ergy E, Eq. (@J is actually a non-linear eigenvalue prob- 
lem. However, to any given order of DPT, E may be 
expanded in a series in J± to obtain an equation with a 
true Hamiltonian form within the degenerate manifold. 
Once \^&o) and E are known, the full wavefunction can be 
reconstructed as |*> = (1 - G^^o) = E,T=o o)- 

To understand the ultimate nature of the ordered state 
on the plateau, we must carry out DPT at least to the 
lowest order at which the degeneracy is lifted. Our ap- 
proach will be to derive the lowest order diagonal (in the 
Sf basis) and off-diagonal terms in J± which remove the 
degeneracy. First, note that any off-diagonal term must 
flip spins in such a way as to preserve the 3:1 constraint 
on each tetrahedron. This can only be accomplished by 
flipping spins around a non-trivial closed loop on the py- 
rochlore lattice. The smallest such loop involves flipping 
spins on 3 different bonds, and flipping a spin from ±| 
to =f| requires Hi to act 3 times, so off-diagonal pro- 
cesses occur first at 9 th order (for general S, they occur 
at 0(J^ S ))! Diagonal processes can occur sooner. The 
first-order splitting vanishes, since the action of Hi cre- 
ates two tetrahedra with spin Sf ^ 3. The energy of the 
plateau state does shift relative to other levels at second 
order. Although it is more non-trivial to see, in fact, 
splitting within the classical manifold does not occur un- 
til sixth order (this is true for any S) ! The reason for this 
is somewhat involved 8] , but follows from the locality of 
Hi and the strong 3:1 constraint. 

By straightforward but somewhat involved calcula- 
tions, one can compute the lowest-order diagonal and 
off-diagonal contributions to the effective Hamiltonian, 
H e fi, explicitly. We find 



H 



eff = a 



ca 9 J Z J2 (Pa) (o B \+ h-c) , (5) 
p 



|*o> = £|*o>, 



(4) 



where the sums are over all hexagonal plaquettes P, 
and the states |Oa), |Ob) represent the two "flip- 
pable" configurations (A/B) with alternating major- 
ity and minority spins around the plaquette. The 
operator £p is a diagonal operator representing the 
6 th order energy splitting of plaquette P. It takes 

the values £ v - fO -2L -9801 s ) ~ 

cue vaiucs c P — (u, 25 , 10Q , 50Q , 625Q ) ~ 

(0,-5.76,-0.27,-1.60,-1.57), respectively for the 

configurations (TTTTTT, 1UUT, ItTtTT, mm, iTTITT) 

(and cyclic permutations), where t / I denote major- 
ity/minority spin states, respectively. This diagonal in- 
teraction is equivalent to a combination of further neigh- 
bor Ising exchange couplings, plus an additional 3-spin 
interaction term; the plaquette representation is, how- 
ever, more convenient and simpler. The constant c = 

53178588 ^10 
12371645 ~ °' 

The Hamiltonian in Eq. j£| (acting within the 3:1 
Hilbert space) is the basis for our discussion of the or- 
dered state on the plateau. Inspection of £p shows that 
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the A/B configurations are significantly more favorable 
than all others. A natural approximation is then to ig- 
nore the energy differences between other configurations, 
which gives 

Hqdm = v5^(|ox)<oa| + |ob>(o b |) 
p 

-Kj2(\OA)(o B \ + h.c), (6) 
p 

where, extrapolating to the physical limit a — » 1, we have 
V « — 5.76J, K — cJ ~ 4.3J, while in the strict small 
a limit, V/K — > — oo. Eq. J^J is useful because it maps 
directly to a type of "Quantum Dimer Model" (QDM)[jfl. 
In particular, the pyrochlore sites can be mapped onto 
the links of a diamond lattice, whose sites lie in the cen- 
ters of pyrochlore tetrahedra. A minority spin can be 
represented as a "dimer" occupying the corresponding 
link, and the 3:1 states become non-overlapping dimer 
coverings of the diamond lattice. The V and K terms 
map directly to standard "potential" and "kinetic" terms 
for these dimer configurations 0. 

QDMs of this form have been considered on a num- 
ber of two and three-dimensional lattices^, Fixing 
K > 0, the ground state depends upon the dimension- 
less parameter v = V/K, and a few general conclusions 
can be drawn. The Rokhsar-Kivelson point, v — 1, is ex- 
actly soluble, and demarks a boundary between ordering 
into configurations with no A/B plaquettes (for v > 1) 
and, for three-dimensional bipartite lattices like the dia- 
mond, a spin liquid phase with no broken symmetry (for 
v a < v < 13], with some critical coupling v c \ which is 
usually positive). For v < v C 2 < i> c i, the ground state is 
again ordered, and adiabatically connected to the diag- 
onal ground state selected by the diagonal V < term 
alone. For the examples (square, honeycomb, and tri- 
angular lattices) in which the phase diagram has been 
studied by quantum Monte Carlo methods [Jjj> the crit- 
ical coupling v C 2 is either slightly negative or even posi- 
tive, v C 2 > —0.2. The persistence of the "diagonal" state 
to small v can be readily understood. For v — > — oo, the 
energy is minimized by the diagonal state which contains 
the maximum density of flippable A/B hexagons. Such 
a configuration is aiso the one most connected by the K 
term to other 3:1 states. Hence, by "order by disorder" 
reasoning, an appropriate superposition "centered" (in 
Hilbert space) around the diagonal ground state above is 
energetically favored by both terms. Explicit variational 
wavefunctions for 7i c ff in this spirit will be described in a 
future long publication. What or how many intervening 
state(s?) might occur for v C 2 < v < v c \ is not known, and 
it is even possible that there might be no intermediate 
state, and instead a direct transition from the diagonal 
ordered state to a spin liquid at v = v c2 — v cl . 

From this reasoning, in the extrapolation to a = 1, 
which gives v w —1.3, we expect the ground state se- 



lected by quantum fluctuations to have the same symme- 
try as that of the easy-axis limit v — > — oo. In particular, 
for the QDM in this limit, we require the 3:1 configura- 
tion with the highest density of A/B hexagons. Because 
these hexagons overlap, and the 3:1 configuration space 
is highly constrained, this is not a trivial problem. It 
is instructive to examine a single primitive unit cell of 
the pyrochlore lattice, which is a region of space enclosed 
by 4 hexagonal plaquettes. By careful inspection, and 
making use of the 3:1 constraint, one finds that at most 
one of these hexagons can be of A/B type. Therefore the 
maximum fraction of A/B hexagons is j. 

By extensive analysis |8|, we have found one single can- 
didate, shown in Fig.Q which saturates this upper bound 
on the density of A/B hexagons, and moreover, gives the 
lowest energy that we have been able to find for the full 
diagonal term (J2p £p) with all plaquette energies in- 
cluded. Having the maximal density of potential Res- 
onating plaquettes, we call it the "R" state. It has a 
magnetic unit cell consisting of 4 structural pyrochlore 
unit cells, and has the symmetries of the P4 3 32 space 
group - particularly significant is that this state pre- 
serves all rotation symmetries of the pyrochlore lattice 
(up to translations). This implies that, contrary to the 
suggestion in Ref. |4|, only isotropic magnetostriction is 
expected for the R state, and not a rhombohcdral distor- 
tion along a (111) axis. We note that a Landau theory 
analysis predicts the symmetry-breaking transition from 
the P4 3 32 space group (R state) to Fd3m (pyrochlore 
symmetry) is first order in agreement with the exper- 
imental conclusions [l^. 

We now turn to other effects of quantum fluctuations. 
Quantum effects lead to a spin gap, A, on the plateau, 
and consequently an activated temperature dependence, 
AF(T) - M z (0) oc exp{-A(h)/k B T). Another quantum 
effect is a suppression of the local "staggered" moment, 
even at T = 0. To leading order in DPT, the differ- 
ence of the moment on majority and minority sites is 
suppressed to (S^^) — (S^ in ) < 2.2, below the classical 

value (S'max) ~~ (^min) = 

Experiments have studied the low-temperature behav- 
ior on exiting the plateau at lower and higher fields, 
h c i , h C 2 , respectively. The lower edge of the plateau for 
both CdCr 2 O4 and HgCr 2 O4 shows a jump in mag- 
netization, which then decreases linearly down to zero 
magnetic field, suggesting the low-field state is adiabat- 
ically connected with the zero field state. The higher 
field edge, observed only in HgCr 2 04, shows a continuous 
transition out of the plateau, and a continuous evolution 
with no further jumps up to the highest fields measured. 
There is however some indication of a "kink" , possibly 
related to the expected lattice symmetry restoration, for 
fields above the plateau but below full saturation. 

The first order/continuous nature of the transitions at 
the lower/upper edges can be explained classically [7], and 
occurs already in the trivial problem of a single classi- 
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cal tetrahedron. Classically, however, above the plateau, 
though the spins become non-collinear, the ordering of 
the transverse moments, if present at all, is critically de- 
pendent upon the nature of further neighbor interactions. 

By contrast, it is determined in the quantum theory 
even for the nearest-neighbor model. In general, a con- 
tinuous quantum transition off the plateau should be de- 
scribed by the "Bose-Einstein Condensation" (BEC) of 
a single magnon excitation with AS Z — +1, above the R 
state, whose gap vanishes at the plateau edge. The result- 
ing state above the plateau has therefore transverse spin 
order, (S^) 7^ 0, in agreement with classical expecta- 
tions. The precise spatial arrangement of the transverse 
ordering is determined by the symmetry of the lowest 
energy magnon, which must form an irreducible repre- 
sentation (irrep) of PA^i2. In the easy axis DPT, the 
magnon is formed from a superposition of states with a 
single minority spin changed from Sf = — § to Sf = — |. 
At 0(a 2 ), the magnon can hop between minority sites, 
and crucially, the effective hopping amplitude is nega- 
tive, despite J > 0, because it occurs by virtual hopping 
through a majority site. The minimum energy magnon 
state then is a uniform plane wave (trivial irrep), which 
implies the transverse spin order is ferromagnetic. That 
is, for h > h C 2 the system has the same space group 
symmetry as the R state. While there is no guarantee 
the trivial irrep remains lowest in energy for a = 1, it is 
our best guess for the nearest-neighbor model. If second 
neighbor exchange J 2 is significant, this conclusion may 
change. In particular, if J 2 is larger than the effective 
hopping amplitude and antiferromagnetic, the lowest en- 
ergy magnon has a different symmetry. The resulting 
transverse order is more complex, with a magnetic unit 
at least 3 times larger than that of the R state 0- In 
any case, the BEC picture implies quite generally that 
the symmetry for h > h C 2 is as low or lower than that 
of the R state. Therefore, since the fully polarized state 
has the full pyrochlore symmetry, with increasing field 
there must be a transition from to Fd3m before 

the saturation field. This may be associated with the 
"kink" observed experimentally. 

The classical model of Penc et aJ.Q successfully ex- 
plains many experimental features of CdCr2 04 and 
HgCr204. However, classically, additional second and 
third neighbor microscopic Heisenberg interactions are 
required to stabilize the R state, as suggested by H. Ueda 
et aI0. Our results show that quantum fluctuations 
generate effective interactions that stabilize the R state 
even when the microscopic exchange is purely nearest 
neighbor. Furthermore, quantum fluctuations may deter- 
mine the symmetry of the transverse spin order expected 
just above the plateau, which, if measured by neutrons, 
would be very telling. Experimentally, the second and 
third neighbor exchanges could be determined by inelas- 
tic neutron measurements of the magnon spectra in the 



fully polarized state, as has been done in Cs9CuCLi|l3|. 
A more direct test of the relevance of quantum fluctua- 
tions would be to look for signatures of BEC criticality 
(as in Refs. fl^fl^h at the upper edge of the plateau in 
HgCr 2 04. We note that H. Ueda et observed an 

increase of h c i with increasing temperature, concluding 
that "thermal fluctuations stabilize collinear spin config- 
urations" . The BEC picture gives an alternative explana- 
tion: the BEC temperature grows as the magnon density 
increases. 

Theoretically, it is intriguing to contemplate the possi- 
bility that the plateau phase might be close to a quantum 
phase transition to a spin liquid state. A field-theoretic 
analysis, based on the mapping of the 3:1 Hubert space 
to a gauge theory with projective symmetry, indeed pre- 
dicts the R state as the simplest phase proximate to the 
spin liquid^. Other calculations, predictions, and com- 
parison to different theoretical approaches will be made 
in a future publication 8] . 
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